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We perform an analytical renormalization group (RG) study to address the role of Coulomb re- 
pulsion, the competition between extended s-wave superconducting order (s±), and the spin density 
wave (SDW) order and the angular dependence of the superconducting gap in multi-pocket models 
of Iron based superconductors. Previous analytic RG studies considered a toy model of one hole 
and one electron pocket. We consider more realistic models of two electron pockets and either 
two or three hole pockets, and also incorporate angular dependence of the interactions. In a toy 
2-pocket model, SDW order always wins over s± order at perfect nesting; s± order only appears 
when doping is finite and RG flow extends long enough to overcome intra-pocket Coulomb repulsion. 
For multi-pocket models, there are two new effects. First, in most cases there exists an attractive 
component of the interaction in s± channel no matter how strong intra-pocket repulsion is, such 
that the system necessary becomes a superconductor once it overcomes the competition from the 
SDW state. Second, in 4-pocket case (but not in 5-pocket case), s± order wins over SDW order even 
for perfect nesting, if RG flow extends long enough, suggesting that SDW order is not a necessary 
pre-condition for the s± order. Our analytic results are in full agreement with recent numerical 
functional RG studies by Thomale et al. arXiv:1002.3599vl] 



I. INTRODUCTION 

Since 2008 a lot of efforts in condensed-matter com- 
munity have been devoted to solve the puzzle of high-T c 
superconductivity (SC) in newly discovered Fe-based su- 
perconductors. To a large extent, in two years the com- 
munity managed to obtain the data for the pnictides in 
the amount comparable to that collected for the cuprates 
over twenty yearsi. 

The family of _Fe-based superconductors is already 
large and keeps growing. It includes doped 1111 systems 
RFeAsO (R =Rare earth element)^ - — , doped 122 sys- 
tems XFe2As2(X=alkaline earth metals )&~—, as well as 
111 and 11 systems like LiFeAaS. and FeTe/Se 10 . The 
parent compounds of most of these materials exhibit a 
spin density wave (SDW) order—, and superconductiv- 
ity emerges upon either hole or electron doping, or upon 
gradual substitution of one pnictide by the other (As by 
P). In some systems, like LiFeAsS- and LaFePO^, SC 
was found already without doping, instead of a magnet- 
ically ordered state. 

ARPESi^, de-Haas van Alphen oscillations measure- 
ments^, and first-principle calculations — ~— all show 
that low-energy electronic structure of pnictides in 2D 
basal plane consists of two nearly circular, non-equivalent 
hole pockets located at the center of the Brillouin 
zone(BZ), and two symmetry-related elliptical electron 
pockets, located near the corners of the BZ in the folded 
zone scheme, or near (0, tt) and (tt, 0) points, respectively, 
in the unfolded zone scheme. [Folded and unfolded zones 
differ in treating the pnictides - folded zone takes into 
account the fact that there are two non-equivalent po- 
sitions of pnictides above and below Fe plane, and has 
two Fe atoms in the unit cell, while unfolded zone incor- 
porates only Fe atoms and has one Fe atom in the unit 
cell.] These hole and electron pockets form warped cylin- 



ders in 3D space. In addition, in some pnictides there is 
a fifth cylindrical hole pocket centered at (tt, tt) in the 
unfolded zone (at (0, 0) in the folded zone, like other two 
hole pockets), while in other pnictides this fifth Fermi 
surface (FS) becomes a 3D sphere centered near k z — tt 
along z— direction. 

A lot of work has been done over the last two years re- 
garding the symmetry of the order parameter and the in- 
terplay between SDW and SC orders. Most of researchers 
(but not all, see Ref. [Lsl ) believe that the gap symmetry 
is extended s— wave (s±), meaning the gap transforms 
according to A\ g representation of tetragonal sym- 
metry group, but the average gap values along hole and 
electron Fermi surfaces have different sign. However, the 
structure of the gap is still a puzzle. Early works based on 
either spin-fluctuation scenario 17 i 19 ~ — or on renormal- 
ization group (RG) study of a toy model of one hole and 
one electron F S 22 i 23 found a simple angle-independent 
s± gap. Subsequent more sophisticated numerical stud- 
ies, which take into account multi-orbital nature of low- 
energy excitations in the pnictides, have reported an- 
gular dependence of the s± gap, with cos2<f> variations 
on the electron FSs and cos4(f> variations along the hole 
FS a 20 i 24 ~ — . The cos2(j) modulations of the s± gap on 
the electron FSs has also been obtained in the analytical 
studj*2£. If cos 2<j) variation is strong enough, the gap has 
"accidental" nodes along electron FSs, still preserving s± 
symmetry. 

Other recent theory proposals include s ++ stated, s± 
state with nodes on hole FSs due to strong cos 40 modu- 
lations^, and s H state with nodes at particular k z along 
z— direction^. 

Experiments are generally consistent with s± gap sym- 
metry, but whether or not the gap has nodes in particular 
materials is still subject of debate^Ir— . In addition, there 
is no information yet from the experiments where the gap 
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nodes are located, if present. ARPES measurements of 
the gap along hole FSs (taken at fixed k z ) on various 
Fe— pnictides^r— indicate that the gap is almost angle 
independent , but the detailed measurements of the gap 
separately along each of the two electron FSs are still 
lacking, with only few exceptions 4 ^. 

From theoretical perspective, the most relevant issue 
is the nature of the pairing interaction. Conventional 
electron-phonon coupling is always a candidate, particu- 
larly when the gap has s— wave symmetry, but has been 
shown to be rather weaki£, and is incapable to account 
for T c ~ 50K even if one neglects destructive effect of 
Coulomb interaction. This leaves electron-electron inter- 
action (i.e., dressed Coulomb repulsion at a finite mo- 
mentum transfer) as the dominant pairing interaction. 
Such interaction cannot give rise to a constant s— wave 
gap, but it can give rise to either momentum-dependent, 
sign reversing gap along a given FS, like d x i^ y i gap in 
the cuprates, or to gaps of different signs along differ- 
ent FSs. Coulomb interaction at large momentum trans- 
fer contributes to the pairing both directly and by cre- 
ating effective pairing interactions mediated by collec- 
tive excitations in either spin or charge channel. The 
close proximity to magnetism makes spin fluctuations a 
preferable candidat e - 17 ' 19 ' 24 , although orbital fluctuations 
were also considered recently^ 7 -. Both direct Coulomb in- 
teraction (the pair hopping) and magnetically-mediated 
interaction are attractive for s± gap, and the total at- 
tractive pairing interaction is a combination of the two. 
Still, to give rise to a pairing, this combined attrac- 
tive pairing interaction has to overcome repulsion coming 
from Coulomb interaction at small momentum transfers. 
A conventional McMillan- Tolmachev renormalization 49 
does not help here because both repulsive and attrac- 
tive parts of the interaction renormalize in the same way, 
and if repulsive part is initially stronger, the renormal- 
ization just reduces the strength of the total repulsive 
interaction, but cannot change its sign. 

How to overcome Coulomb interaction became the ma- 
jor issue for s— wave pairing in the pnictides. The situa- 
tion is further complicated by the fact that low-energy ex- 
citations are composed of all 5 hybridized Fe 3d-orbitals, 
and Coulomb interactions between fermions belonging to 
a given orbital and belonging to different orbitals are 
equally important. Only if all of these interactions are 
approximated by the same momentum-independent Hub- 
bard U, this U cancels out in the pairing problem and 
one is left with a pure spin-mediated interaction (this was 
termed "Coulomb avoidance"- ). In reality, Coulomb in- 
teraction is momentum-dependent and is larger at a small 
momentum transfer than at a large momentum transfer, 
and intra-orbital and inter-orbital interactions are also 
different. As a result, direct Coulomb pairing interaction 
is repulsive for superconductivity with angle-independent 
plus-minus gap. Spin-mediated interaction is attractive 
and can potentially compete with direct Coulomb repul- 
sion. However, at least at weak/moderate coupling direct 
Coulomb repulsion is the largest term. This holds even 



when magnetic correlation length £ diverges because for 
the pairing one generally needs the interaction at non- 
zero frequencies, where it remains finite. 

There are two possibilities to obtain s± superconduc- 
tivity despite strong Coulomb repulsion. First, multi- 
orbital character of excitations in the pnictides implies 
that the attractive pairing interaction (either pair hop- 
ping or spin-fluctuation exchange) is angle- dependent. 
The angle-dependence comes from coherent factors which 
dress up the interactions when one transforms from the 
orbital picture -in which different parts of the Fermi sur- 
face are made of different orbitals- to the band picture, 
in which free-fermion part is simply e^ctc^, and all infor- 
mation about multi-orbital character is passed onto in- 
teractions. Once interaction is angle- dependent, the gap 
also becomes angle-dependent, and the system adjusts 
the angle-dependence of the gap to minimize the effect 
of Coulomb repulsion (we discuss this in more detail be- 
low) . The most natural is the case when the gap acquires 
± cos2</> components along the two electron FSs, and the 
magnitude of this component is adjusted to balance the 
interplay between small q Coulomb repulsion and the 
combined attractive interaction in s± channel. When 
Coulomb repulsion dominates, the angle-dependent part 
is large, and the gap has four nodes along each of the 
electron FSs. 

Second, one can analyze how the interactions evolve 
as the system flows towards smaller energies, relevant 
to superconductivity. This flow involves renormal- 
izations of interactions in both particle- hole (p-h) and 
particle-particle (p-p) channel, and goes beyond RPA. 
This flow has been studied numerically, using func- 
tional renormalization group (fRG) technique 2 ^— and 
analyticall y 22 ' 28 ' 51 , within parquet RG. Both are weak- 
coupling studies, based on the Hamiltonian which con- 
tains screened Coulomb interaction, but no additional 
spin-fluctuation interaction. The results of both types 
of studies are quite similar: it turns out that Coulomb 
repulsion at small momentum transfers decreases upon 
system flow to smaller energies, but the pairing interac- 
tion at large momentum transfers (the pair-hopping from 
hole to electron FSs), which is attractive for s" 1 super- 
conductivity increases. The increase of the pair hopping 
is the result of the "push" from the inter-pocket density- 
density interaction which by itself leads to SDW insta- 
bility. If RG flow of the couplings persists over a wide 
enough range of energies, pair-hopping interaction ex- 
ceeds Coulomb repulsion, and the system develops an at- 
traction in the s± channel. In this situation, the ±cos24> 
variations of the gap on electron FSs induced by angle 
dependence of the interaction, are not crucial for the 
pairing, and the system can develop an s± gap without 
nodes. 

While this scenario is quite generic, the earlier parquet 
RG stud y 22 ' 51 was more limited in scope. It was done for 
a toy model of one hole and one electron FS centered at 
(0,0) and {n,Tr), respectively. For such a model, angle 
dependencies of the interactions must be symmetric with 
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respect to interchanging x and y momentum components 
both near k = (0,0) and k — (tt,tt) and can only be in 
the form cos 40, cos 80, etc which are subleading terms 
in the expansion of Ai g gap for a single FS (no cos 20 
terms!). Such cos4n0 terms are generally irrelevant and 
were neglected in toy model analysis, i,e. the gap along 
each of the FSs was approximated by a constant. Like 
we said, for momentum-independent gaps, the bare in- 
teraction in the s H channel is repulsive if intra-pocket 
Coulomb repulsion is the largest. The interaction flows 
under RG and changes sign at some value of RG pa- 
rameter. Still, all along parquet RG flow, the pairing 
interaction remains secondary to the interaction in the 
SDW channel. As a result, for perfect nesting the sys- 
tem develops an SDW order. Only when the system is 
doped and the logarithmical flow of the SDW vertex is 
cut at low energies, SC channel takes over and the system 
develops an s H superconductivity. 

In the present paper we extend earlier parquet RG 
analysis to multi-pocket models of Fe-pnictides. We con- 
sider two models. The first one has two electron FSs, 
located at (0, tt) and (tt, 0) in the unfolded zone, and two 
hole FSs located near (0,0). The second model has an 
additional hole pocket centered at (tt, tt) in the folded 
zone. The presence or absence of this additional FS in 
different Fe-pnictides is attribute d 19 ! 27 to the difference 
in the distance between the pnictide (e.g., As or P) and 
the Fe-plane. To avoid overly complicated analysis of RG 
equations we only consider two limiting cases: one when 
the two hole pockets centered at (0, 0) are completely 
equivalent, and the other when one pocket is coupled to 
electron FSs much weaker than the other one and can 
be neglected [this effectively reduces 4 pocket model to 
3 pockets and 5-pocket model to 4 pockets]. The sys- 
tem behavior is identical in the two limits which gives us 
confidence that it remains the same also in between the 
limits. Throughout the paper we assume that all FSs are 
cylindrical and neglect their variations along k z . 

We argue that new physics: nodal s ± SC, appearance 
of SC even at perfect nesting, emerges once one extends 
the model from 2 to 4 pocket (or to 3-pocket when only 
one hole FS is relevant). This result is consistent with 
early assertion a 22 ' 28 that the 3-pocket model is the min- 
imum model needed to understand the key physics of 
Fe— pnictides. The behavior of 5 pocket model, on the 
contrary, is in many respects similar to that for 2-pocket 
model (except that nodal s H SC is still possible). That 
4-pocket and 5-pocket models behave differently under 
RG has been the conclusion of fRG study by Thomalc et 
dl (Ref.,27), and our results fully agree with their analysis. 

We extend previous 2-pocket study in two directions. 
First, we incorporate cos 20 angular dependence of the in- 
teractions which give rise to ± cos 20 modulations of the 
gaps on electron FSs. In this situation, there are at least 
three different effective vertices for s± gap symmetry, and 
we argue that, if the dominant angle dependence comes 
from electron-hole interaction, one of them remains pos- 
itive (i.e., attractive) over the entire RG flow even when 



bare intra-pocket repulsion is the largest interaction. We 
show that the stronger is the angular dependence of the 
interaction, the stronger is the tendency to develop an 
nodal s± order. 

Second, we re-analyze the interplay between SC and 
SDW channels. For angular-dependent interactions, are 
also several SDW vertices of which at least one is attrac- 
tive along the whole RG trajectory. We compare the flow 
of the leading vertices in SDW and SC channel. We show 
that in 4-pocket model, the trajectory of the leading SC 
vertex are steeper than that of the leading SDW ver- 
tex, and the latter remains larger only down to some RG 
scale. At smaller scales (i.e., when RG flow extends fur- 
ther to lower energies) the SC vertex overshadows SDW 
vertex even at perfect nesting. This agrees with fRG 
study by Thomale et .al^. We argue, based on our an- 
alytic consideration, that the crossing between SC and 
SDW vertices under RG flow is to a large extent a com- 
binatoric effect - compared to 2-pocket case (where SDW 
and SC vertices flow to the same value under RG), the 
presence of the two electron FSs adds the factor of 2 to 
the renormalization of the SC vertex as the pair hopping 
can, e.g., hop a pair of k, —k fermions from the hole FS 
to each of the two electron FSs, but momentum conser- 
vation does not allow such factor of 2 to appear in the 
renormalization of the SDW vertex. 

We further find that for 5-pocket model, SC vertex 
always remains secondary compared to SDW vertex, just 
like in 2-pocket model. Furthermore, like in 2-pocket 
model, SC and SDW vertices flow to the same value at 
the fixed point. This again agrees with fRG result by 
Thomale et. at. 

The rest of the paper is organized as follows: Section 
HTl outlines our approach, explains the subtleties in the 
RG flow, and discusses the technique by which we in- 
corporate the momentum dependence into our analysis. 
In Sec. IIIII we briefly review the results for the 2-pockct 
case. Section HVl is the central section in the paper - here 
we consider in detail the 4-pocket model. We discuss 
SC vertices in the presence of angular dependence of the 
interaction, the RG flow, and the competition between 
SC and SDW instabilities. Most of our treatment in this 
section is for the limit when one of the hole FS can be 
neglected. Later in the section we analyze another limit 
when the two hole FS are equivalent and show that the 
system behavior is identical in the two limits. In Sec. |V] 
we discuss 5-pocket model and argue that its behavior 
to a large extent is similar to that in 2-pocket model. 
We summarize our results in the Conclusion. 



II. DISCUSSIONS 

Before starting the detailed description of each model, 
we present a brief discussion on the outline of our ap- 
proach. We first describe the central idea behind our 
analysis, then describe the technique through which we 
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incorporate the angular dependence of the vertices and 
finally discuss how we incorporate these into RG equa- 
tions. 



A. Approach 



We consider four-fcrmion interactions between 
fermions located close to the FSs of two or more pockets. 
We consider Hamiltonian with all possible quartic 
interactions allowed by symmetry and ask what can 
be said about the onset of SC, SDW, and, possibly, 
CDW instabilities. The usual approach is to write 
down equations for effective vertices Ti in SDW, SC, 
CDW channels and check for the existence of critical 
temperatures T^ s at which I\ diverges. In case of 
competing instabilities, the equations for the effective 
vertices are coupled and, once the coupled system is 
solved, the instability with the highest Ti ns {T c ) takes 
over. 

As we said in the Introduction, when one does such an 
analysis, one deals with interactions taken at the scale of 
T c , which are not the same as the terms in the Hamilto- 
nian. To account for the flow of the couplings from the 
scale of the bandwidth down to T c , we need RG analy- 
sis. This analysis assumes renormalizability of the the- 
ory and can be rigorously justified only when the RG 
flow is logarithmical (i.e., interactions vary as functions 
of the logarithm of the running scale E). One well-known 
example of logarithmical RG flow is the renormaliza- 
tion in the particle-particle channel (Cooper renormal- 
ization). Another, specific to our case, is the renormal- 
ization in the particle-hole channel, involving intermedi- 
ate fermions from hole and electron pockets. Because 
hole and electron dispersions are of opposite sign, such a 
rcnormalization also generates logarithmical dependence 
of the running energy and/or momentum as long as the 
running energy exceeds the sum of energies of the top of 
the hole band and the bottom of the electron band (hole 
and electron masses do not have to be equal). 

The logarithmical renormalizations in the particle- 
particle and particle-hole channels are characterized by 
corresponding polarization bubbles. Let c describe a 
hole band centered at k — and / describe an electron 
band centered at Q. Assume for simplicity that hole and 
electron masses are equal. For a perfect nesting, hole 
and electron dispersions obey e c (fc) = eo — fc 2 /(2m/,.) = 
—Ef(k + Q). The two logarithmically singular polariza- 
tion bubbles are 
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where the dots stand for non-logarithmic terms, E = 
max{£l, vfq} and E > Ep, A is the upper cutoff of or- 
der bandwidth, and the propagators are given by G x = 
- — — =■, x being c or f. 

The RG study requires caution as the couplings flow 
differently for energy scales above Ep and below Ep. The 
reasoning is simple: logarithmical RG analysis requires 
that internal momenta in each diagram for vertex renor- 
malization be larger than external momenta, which are 
of order kp. When typical internal energies are larger 
than Ep, internal momenta are larger than kp, and ver- 
tex corrections in both particle-particle and particle-hole 
channel are logarithmic. This gives rise to parquet RG. 
When typical energies are smaller than Ep, the strength 
of the renormalization in the particular channel depends 
on the interplay between external momenta. When to- 
tal incoming momenta is zero, renormalization in the 
particle-particle channel is still logarithmically singular, 
but in the renormalization of the particle-hole channel, 
the logarithm is cut by external Ep. Conversely, for 
the vertex with transferred momentum equal to the dis- 
tance between hole and electron FSs, the renormalization 
in the particle-hole channel is still logarithmical, but in 
the renormalization in the particle-particle channel, the 
logarithm is now cut by external Ep. As a result, the 
renormalizations in the particle-particle and particle-hole 
channels are coupled at energies above Ep, but become 
decoupled at energies below Ep. At E < Ep parquet RG 
equations are replaced by conventional ladder RG equa- 
tions dTi/dl = Ff, where I = log Ep/E. Thus the flow 
of the couplings splits into the flow from the bandwidth 
down to Ep, where different vertices are all coupled, and 
the flow below Ep, where different vertices are decoupled 
(see Fig. [1]). This reasoning is particularly important in 
our case, as in pnictides Ep ~ lOOmeV is much smaller 
than the bandwidth, which is a few electron volts. 

Depending on the character of the flow, the bare values 
of the couplings, and the ratio of the bandwidth and the 
Fermi energy, several situations are possible and shown 
in Fig. U 

• The RG flow diverges and the normal state becomes 
unstable before the scale of Ep is reached (Fig. 
. In this situation, the instability is reached 
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FIG. 1: Illustration showing how the couplings evolve under 
RG flow. The bare couplings (the parameters of the Hamil- 
tonian) are defined at energies comparable to the bandwidth 
A. In pnictides, this scale is 2 — 3eV, much larger than the 
Fermi energy Ef ~ O.leV SC and SDW instabilities likely 
come from even smaller energies because instability temper- 
atures are at least order of magnitude smaller than Ef- The 
couplings vary as one integrates out higher energies. This 
variation (i.e., the flow of the couplings from A down to the 
running scale E) can be generally captured by the RG tech- 
nique. In pnictides, the flow is different above and below Ef- 
Above Ef, each of the couplings changes because of inte- 
grating out higher energy fermions in both particle-hole and 
particle-particle channels. The RG equations in this region 
are called parquet RG, because renormalizations extend in 
the two directions. Below Ef, each vertex continue to flow 
due to renormalizations in only one channel, either particle- 
hole or particle-particle, depending on the external momenta. 
The RG equations in this region is called ladder RG, because 
renormalizations extend only in one direction. 



already within parquet RG. This situation is the 
most interesting one from physics perspective par- 
ticularly because in several cases different vertices 
diverge simultaneously, and the fixed point has an 
enhanced symmetry (e.g., in the 2-band model, the 
vertices in SDW, SC, and orbital CDW channels 
all diverge at the fixed point which then has 0(6) 
symmetr y 51 ' 52 ). The instability at energies above 
Ep is, however, unlikely scenario for the pnictides 
because the largest instability temperature is only 
a fraction of Ep. 

• The RG flow reaches Ep before couplings diverge 
(Fig. [5Jd). In this situation, parquet RG creates a 
hierarchy of the couplings at Ep: Ti(Ep). Below 
Ep, different Ti decouple and, for a perfect nest- 
ing, each continue evolving according to a ladder 
RG, i.e., like T^E) = T^Ep)/^ - T^Ep) log ^f). 
The instability occurs at the energy (temperature) 
at which Ti(Ep) log ^ = 1. Obviously, the win- 
ning channel is the one in which the coupling is the 
largest at Ep. 

• When nesting is not perfect (i.e., energies of the 
top of the hole band and of the bottom of the elec- 
tron band are not exactly opposite), the coupling 
in the SC channel continues to follow T SC (E) = 
T sc (E F )/{l-r sc (EF) log ^) simply because SC 
instability involves pairs of fermions with k and — k 
from the same FS and is non-sensitive to a devia- 



tion from perfect nesting. However the logarithmic 
flow of SDW and CDW vertices is now cut at some 
scale Eb. Suppose Eb < Ef (Fig. [2k). In this sit- 
uation, SC eventually wins over SDW and CDW 
instabilities, even if superconducting T is sublead- 
ing at Ep (but it needs to be attractive at Ep). 

• When Eb exceeds Ep (Fig. [2ji) particle-hole and 
particle-particle channels decouple already within 
the applicability range of parquet RG. At E < 
Eb, SC vertex continue to grow as T SC '(E) = 
T sc {E b )/{l - T sc (E b ) log E b /E), while vertices in 
density-wave channels get frozen at their values at 
E ~ Eb. In this situation, SC instability again 
wins, even if it was subleading at Eb, provided that 
the superconducting vertex is attractive at Eb. 

A subtle point: below we will be presenting the RG 
flows of vertices and couplings at energies both above and 
below Ep in terms of the logarithmic variable L ~ logE. 
One has to bear in mind, however, that the prefactor 
for the logarithm actually changes between E > Ep 
and E < Ep because for E > Ep the integration 
over intermediate energies involves only positive excita- 
tions for electron states (and negative for hole state), 
i.e., j — (to/27t) Jp de k , while for E < Ep, one 

has to linearize the dispersions of holes and electrons 
near the FS and integrate on both sides of Ep, i.e., 

I 



d z k 



(m/n) f^! F de k . To simplify the presentation 



(2/r) 2 

we just define 



L = 



Un± ,E>E F 
Un^h+lnSf ,E < E F 



and use the same symbol L for all energies. This is valid 
as long as we describe the system behavior at E Ep 
and E <C Ep. The behavior in the crossover regime 
E ~ Ep is more complex, but this is beyond the scope 
of the present paper. We will also use Lp F = \ln-^. 



B. Incorporating angular dependence 

As we said in the Introduction, multi-orbital character 
of low-energy excitations in Fe-pnictides implies that in- 
teractions between fermions located near hole or electron 
FSs depend on the angles along the FSs. To obtain angu- 
lar dependence of various couplings from first principles, 
one has to transform from five-orbital to five-band picture 
and dress up Coulomb interactions by coherence factors. 
This has been done in several studies (see, e.g., Refs. 
I24II26II53I ) , under the assumption that the Coulomb inter- 
action is so strongly screened that it can be replaced by 
Hubbard U. This assumption is generally valid for sys- 
tems with large FSs, because of many available particle- 
hole pairs for screening, but in systems with small FSs, 
like pnictides, much fewer number of particle-hole pairs 
are available, and the screening of Coulomb interaction 
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FIG. 2: Schematics of the RG flow of SC and SDW vertices T SDW and T 3± . The horizontal scale is L = \ln± for E > E F 
and \ln£- + In^g- for E < E F . SC and SDW vertices remain coupled at energies larger than Ef but decouple below Ef- 
Depending on the bare values of the vertices, the ratio A/Ep, and the doping, four different scenarios are possible (panels 
(a)-(d)). For perfect nesting, there are two possibilities: (a) the vertices diverge at the same scale before the scale of Ef is 
reached. The ratio of the vertices not necessarily tends to one, though, (b) Ef is reached before the vertices diverge. Then, 
below Ef, the vertices decouple, flow and diverge independently, each on its own scale. The vertex that had larger value at 
Ef diverges first and sets the instability. For non-perfect nesting (e.g., at a finite doping), SDW vertex eventually does not 
diverge. SC vertex still diverges, and the system becomes a SC even if SC instability was subleading at perfect nesting. The 
flow of the SDW vertex levels off either at Et < Ef (panel c), or at Et > Ef (panel d). 



is much weaker, and is actually a rather non-trivial phe- 
nomenon at small kp^ Because of this complication, the 
"first-principle" analysis of the angular dependence of the 
interaction is a rather difficult task. 

One can, however, attempt to extract these angular de- 
pendence from symmetry considerations, like it has been 
done for the cuprates^ This is what we will do. Consider 
first the pairing vertex between fermions with incoming 
momenta k and — k and outgoing momenta p and —p. 
Quite generally, for tetragonal D^h symmetry group, this 
interaction can be divided into one-and two-dimensional 
representation, and one-dimensional representation can 
be further divided into Ai g , B\ gi i?2g, and Ai g harmon- 
ics, depending on the symmetry under the transforma- 
tions under k Xjy — > —k x ,y and k x — > k y . Basic func- 
tions from different representations do not mix, but each 
contains infinite number of components, s— wave pair- 
ing corresponds to fully symmetric A\ g representation. 
The s— wave pairing interaction can be quite generally 
expressed as 

U(k,p) = J2 Am n^m{k)^n(p) (2) 

m,n 

where \I/ m (fc) are the basis functions of the Ai ff sym- 
metry group: 1, cosk x cosk y , cosk x +cosk y , etc, and A mn 
are coefficients. Suppose for definiteness that k belongs 
to hole FS and is close to k — 0. Expanding any wave 
function with Ai g symmetry near k = 0, one obtains 
along |k| = kp, 

*m(fc) = a m + b m cos40 fe + c m cos80 fc + ... (3) 

If p is near the same hole FS, the expansion of 
* n (p) also involves cos4(/) p , cos8</> p , etc. There are 
no fundamental reasons to expect that b m , c m , etc are 
much smaller than a m , but sub-leading terms are of- 
ten small numerically. Two known examples are the 



numerical smallness of cos6</), etc components of the 
d x 2_ y 2— wave gap for spin-fluctuation mediated pairing in 
the cuprates^ and the numerical smallness of cos 4<fi, etc 
components of the gap along the hole FSs in fRG— and 
RP A 24 ' 61 calculations for 5-band Hubbard-type model for 
the pnictides. Taking these examples as circumstantial 
evidence, we assume that cos4</>, etc terms are small. If 
so, the interaction between fermions belonging to the hole 
FS can be approximated by angle-independent term. 

The situation changes, however, when we consider 
pairing interaction between fermions belonging to dif- 
ferent FSs. Suppose that k are still near the center of 
the Brillouin zone, but p are near one of the electron 
FSs, say the one centered at (0,7r). Consider all possi- 
ble ^ n (p) with Ai g symmetry. A simple experimentation 
with trigonometry shows that there are two different sub- 
sets of basic functions: 

A:l, cosp x cospy, cos2p x + cos2p y ... 

A : cosp x + cospj,, cos3pa; + cos3pj,... (4) 

Functions from class A have the same properties as be- 
fore - they can be expanded in series of cos Al(f> p (I is inte- 
ger). Functions from class A are different - they all van- 
ish at (0, tt) and are expanded in series of cos(2^> p + 4Z0 p ) 
(i.e., the first term is cos2<^ p , the second is cos60 p , etc). 
Let's make the same approximation as before and neglect 
all terms with I > 0. The functions from class A can then 
be approximated by a constant, but the functions from 
class A are approximated by cos 2(f> p . As a result, s— wave 
pairing interaction involving fermions from hole and one 
of the two electron FSs (labeled as e\) has a generic form 
of 



u euh (k,p) = U eu h +ttei,/iCOs20 Pel + .... 

= u ei)fe (l + 2acos2^ Pel ) + ... (5) 
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where dots stand for cos40fc, cos40 p , cos60 p , etc 
terms. We emphasize that the constant term and the 
cos 2cf>p term in ([5]) are the leading terms of the two sub- 
sets of interaction terms, each form series in cos40fc jP . 
By the same reasoning, the interaction between fermions 
near two electron FSs centered at (0, it) and (n, 0) is ex- 
pressed as 

u ei ^ e2 (k,p) — u ei ,e 2 (1 + 2a' (cos20fe cl + cos20 Pe2 ) 

+4a cos20fc el cos20 Pc2 + ... (6) 

Observe also that the cos 20 terms in ([5]) and © 
change sign under the transformation x — > y (like d x 2_ y 2 
interaction in the cuprates), hence the prefactor for 
cos20 p term in ([5]) changes sign between the two elec- 
tron FSs [cos20 Pel — > — cos20 Pc2 ]. 

The pairing interaction in the form of Eq. ([5]) has 
been introduced in Ref. H|| The authors of Ref. I2H 
however, didn't include into consideration the fact that 
band description is obtained from multi-orbital descrip- 
tion and argued that a must generally scale as k 2 F and 
should be small when kp is small. In fact, the angular 
dependence produced by the coherent factors associated 
with the hybridization of 5 Fe bands are not small in kp 
(Refs. H4II53I), hence a, a', a do not have to be small. 
Accordingly, we will keep a's as just parameters. 

Once the pairing interaction has the form of Eqs. ([5]) 
and ©, the gap along hole FS is still angle- independent, 
but the gaps along the two electron FSs are in the form 
A e ± A e cos 20. When A e is small compared to A e , the 
gaps on electron FSs are nearly angle-independent, but 
when |A e | > |A e |, they have nodes at "accidental" values 
of 0. 

C. The RG analysis with angle-dependent 
interactions 

The angular dependence of the interaction is the key 
element of fRG approach, and this approach uses the full 
momentum dependence of the interactions, i.e., the full 
series of cos 4/0 and cos(20 + 4/0) terms. At the same 



time, fRG approach assumes renormalizability (i.e., that 
the right-hand side of each fRG equation contains only 
renormalized couplings, not some combinations of bare 
and renormalized couplings) . In our analytical approach, 
we do calculations within logarithmic approximation in 
which case we explicitly preserve renormalizability. 

We found, after explicitly evaluating the renormaliza- 
tions of angle-dependent vertices, that the only way to 
justify RG procedure in this situation is to keep the an- 
gular dependence away from RG flow, i.e. allow u ei:C21 
u e i.h, and u e 2,h to flow under RG, while a, a' and a are 
kept unchanged. This can be rigorously justified when 
a's are small and all terms of order a 2 are neglected. We 
will assume without proof that the results of RG analysis 
are valid even when a < 1. There are no new physical 
effects at a < 1 compared to a << 1, so the results at 
a < I should be at least qualitatively correct by conti- 
nuity. 

As we said in the Introduction, the main goal of our 
analysis is to understand whether there are qualitative 
differences between RG flow and the pairing in 2, 4, 
and 5 pocket models for Fe— pnictides. The next three 
sections deal with the comparison of 2,4, and 5-pocket 
models. 



III. THE 2 POCKET MODEL 

This has been studied before ( Refs. I22ll5il ) and we 
briefly review it here to set notations and for further com- 
parisons with multi- pocket cases. 

The 2 pocket model is a toy model consisting of one 
hole pocket in the center the folded Brillouin Zone (BZ), 
and four electron pockets at the corners as shown in Fig 
El To obtain parquet RG equations, we consider ener- 
gies larger than Ep and E^. At such energies deviations 
from perfect nesting become irrelevant, and we can set 
Ep, Eb — > and take hole and electron dispersions to be 
just opposite in sign. 

The interaction Hamiltonian for the 2-pocket model is 
(following earlier notations 2 ^,) 



2tt 



-Hi, 



ui c iisfi 2 s'fp 



y M 2 C pis/jI 2 s' C P4s'/p 



EU3 
2 



V pi 



fp4S'fp 



h.c. 



(7) 



/ j r > Jp 1 sJp2s'JPiS'Jp3S^2^, O P1S P2S' P* S C P3S 



where m are dimensionless couplings, and m is quasipar- 
ticle mass (e f(k) = k 2 /2m — [i). The sum is over the spin 
indices s and s' and the vector momenta pi,p2,P3,P4 with 
momentum conservation assumed. The c and / fermions 
reside at the hole and the electron bands respectively. We 



remind the reader that there is no angular dependence of 
the interactions here because the first angular term that 
comes in is cos 40 which we ignore in our approximation. 
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FIG. 3: Hole (center) and electron (corners) FSs in the folded 
BZ for 2-pocket model. The arrows with symbols indicate 
intra-pocket and inter-pocket pairing interactions (114 and 115 
are intra-pocket interactions, and M3 is inter-pocket interac- 
tion). There also exist density-density and exchange interac- 
tions between hole and electron pockets ( Hi and U2 terms, 
respectively, not shown). 



A. The Vertices 

We begin by looking into the vertices in SC and SDW 
channels. For this, we introduce infinitesimally small SC 
and SDW order parameters A° sc and A° SDW , dress them 
up by including multiple interactions as shown diagram- 
matically in Fig. gj and write the renormalized order 
parameters in the form 



obtain from Fig. gj T SDW 
we obtain from Fig. @] 



U1+U3. For the SC channel, 



1 — U5L —U3L 
-u 3 L 1 - u 4 L j \ A° 



A, 



(9) 



where A e are order parameters on hole and electron 
FSs. Diagonalizing this set and casting the result in the 
form of Eq. |(5J), we obtain two SC Vs. One corresponds 
to a conventional s— wave pairing, is repulsive for all pos- 
itive Ui and is of no interest to us, another corresponds 
to s± pairing and is given by 



For 114 = us, T s± reduces to T 



s± _ 



-Ui + U3. 



The SDW vertex is attractive for positive U\ and u 3 , 
while r s± is attractive only when inter-band pair hopping 
term exceeds intra-band repulsive interaction. Like we 
said, this is very unlikely because both interactions orig- 
inate from Coulomb interaction, and screened Coulomb 
interaction at small momentum transfer (i.e., U4 and u$) 
is larger than that at large momentum transfer (i.e., 1*3). 

To understand whether the negative sign of T s± can be 
reversed, we need to consider RG flow of the couplings. 
This what we do next. 



A i = A?(l+IVL) (8) 
where Ti satisfies ^ = T? 

For a given i, Aj becomes nonzero even for vanishing 
A° when the corresponding Tj diverges. 




FIG. 4: Diagrams for the renormalization of infinitesimally 
small SDW and SC vertices, added to the Hamiltonian in 
order to calculate response functions. Unshaded triangles - 
bare vertices, shaded triangles - full vertices, wave lines - 
fully renormalized interactions. Solid lines correspond to c 
fermions and the dashed line to / fermions. Ai is SDW vertex 
and Ah and A e are SC vertices on hole and electron FSs. 

The computations for the 2-pocket model is straight- 
forward. For the SDW order parameter, we immediately 



B. RG flow between A and Ef 

The RG equations for the couplings have been obtained 
ir*22, and we just quote the result: 

U4 = — [u\ + u 3 ] 

u 5 = - [u% + uj] 

111 = \ux + U3] 

112 = [2uiU2 — 2^2] 

u 3 = [4uiu 3 - 2u 2 u 3 - [u4 + u 5 )u 3 ] (11) 

The derivatives are with respect to L. These equations 
have a single non-trivial fixed point at which all cou- 
plings diverge and tend to u 3 = \/5ui,U4 = u§ — —u%, 
U2 oc (tii) 1 / 3 . The flow of SDW and SC vertices is shown 
in Fig. |6ji. In the process of RG flow, the SC ver- 
tex r s± changes sign and become attractive. The ra- 
tio Y s± /T SDW remains smaller than one during the flow, 
but tends to one upon approaching the fixed point, i.e., 
if this fixed point is reached within parquet RG, super- 
conducting and SDW instabilities occur simultaneously, 
and the system actually cannot distinguish between the 
two. There is another vertex which tends to the same 
value as SDW and SC vertices - it corresponds to an 
CDW instability with imaginary order parameter (an in- 
stability towards orbital currents), The combination of 
3-component SDW, 2-component SC and 1-component 
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CDW instabilities makes the fixed point 0(6) symmet- 
ric^. 

The sign change of the superconducting r s± is the 
most notable effect within the parquet RG flow. Its 
physics originates in the effective "attraction" between 
SC and SDW fluctuations (not the order parameters!), 
namely from the fact that U3, which is the attractive 
component of T s± , gets the boost from m, which con- 
tributes to Y SDW . The boost is 4uiU3 term in the r.h.s. 
of the RG equation for ii 3 . This term overshadows the 
negative effect from it 2 , W4, and u 5 , and as a result 113 
increases under RG. At the same time, intra-pocket re- 
pulsions 1*4 and U5 decrease under RG. At some point, 
u\ becomes lager compared to U4M5, and T s± becomes 
positive. 

However, as we said earlier, parquet RG is only valid 
at energies above Ep. It is unlikely that the fixed point 
is reached above Ep, otherwise there would be at least 
pseudogap effects present above Ep, but there is no 
strong evidence for pseudogap in the pnictides. More 
likely, the couplings evolve under parquet RG (and r s± 
possibly changes sign at some scale) , but Ui remain finite 
at Ep. To continue below this scale we need to derive 
a different set of equations, for which Ui(Ep) serve as 
initial conditions. 



(0) 

Pi 



(b) 



p 2 
p" 

k 



P 3 +Q 



P 4 +q 



(a) 



k+Q 
P+Q 



P ] 


, — » — 

p+Q 


k J 
—§---. 

— 


* k+Q 

5 s ' 




p : 


> k+Q 


f J 


\ -k+Q 

* * — 



FIG. 5: (0) The«3 vertex with general momenta pi, p2, P3+Q, 
P4 + Q (all pi are small and P1+P2 = P3 +Pi)- During calcu- 
lations, three kinds of U3 vertices arise-(a) the one with pi = 
P3, (b) the one with p2 = pz, and (c) the one with pi +P2 = 0. 
The vertex 'b' contributes to the renormalization in p-h chan- 
nel, and the vertex 'c' contributes to the renormalization in 
the p-p channel. 



C. RG flow below the scale of Ef 

The RG equations below Ep for 2-pocket model have 
been derived in Ref. [5l| and we just quote the result. 
The most essential difference with the previous subsec- 
tion concerns 1/3 vertex, which contributes to both SC 
and SDW channels. Below Ep the structure of the ex- 
ternal momenta becomes relevant, and one has to dis- 
tinguish between U3 ' with zero transferred momentum, 

4 with momentum transfer Q, and 4 with zero to- 
tal momentum (see Fig. Fig. [3]). Each of the vertices 



now undergoes logarithmic renormalization in its own 
channel, crossed renormalizations no longer contribute 
because internal E — O(Ep), and the arguments of the 
corresponding logarithms become 0(1). The new equa- 
tions are 



(a) 



.(a) 
M 3 


= 2iti4 — 2uiu 


■(b) 
«3 


= 2iti4 


• (c) 
u 3 


— — [«4 + U 5 ]u^ 


«4 


= -[ul + {uff] 


115 


= -[^+(4 C) ) 2 ] 


iii 




u 2 


= [2uiit2 — 2^2] 



(12) 



where the derivatives are with respect to L. The effective 
vertices in the SDW and SC channels also get modified 
and become 



V SDW = Ul+u W 



> 4 + «B) + \J(U4 -U 5 ) 2 +4(4 C) ) 5 



One can easily verify that new vertices satisfy 



dL 



(13) 



(14) 



as it should be because SC and SDW channels are now 
decoupled (no cross-terms in RG equations). 

The behavior of the vertices below Ep is illustrated in 
Fig. [B]b,c. If SC vertex is already positive (attractive) at 
Ep (Fig. [6p), it diverges at some scale below Ep, but for 
perfect nesting SDW vertex diverges first. Upon doping, 
SDW vertex levels off, and the first instability eventually 
becomes the SC one. If SC vertex remains negative at 
Ep (Fig. [6t), it decreases below Ep but still remains 
negative. In this situation, SC does not appear even 
when SDW order gets killed by non-nesting. 

To summarize: in a 2-pocket model three scenarios 
are possible: (i) the instability occurs simultaneously in 
SDW, SC, and CDW channels, and a fixed point has 0(6) 
symmetry, (ii) SDW instability wins at prefect nesting, 
but yields to s± superconductivity upon doping, and (iii) 
SDW instability exists near perfect nesting, but no SC 
instability emerges when SDW order is suppressed by 
doping. For the cases (i) and (ii), the SC gap has a simple 
plus-minus form, i.e., the gaps along hole and electron 
FSs are angle independent (up to cos 4</> terms which we 
neglected) and are of the opposite signs. 

The 2-pocket model is indeed only a toy model for the 
pnictides. The actual band structure of the pnictides in- 
cludes two electron FSs and at least two hole FSs. The 
question we now address is whether qualitatively new be- 
havior emerges when we increase the number of pockets. 
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We argue below that there are new features not present 
in a 2- pocket model. 

IV. 4 POCKET MODEL 

We now consider the case of 2 hole and 2 electron FSs. 
We neglect k z variation of the FSs and consider a cross 
section in XY plane. The two electron FSs are generally 
ellipses, centered at (0, 7r) and (tt, 0) in the unfolded zone. 
The two hole FSs are circles centered at (0,0). They 
generally are of non-equal sizes, and one is less nested 
with electron FSs than the other. 

The RG analysis of a generic 4-pocket model is 
straightforward but rather cumbersome. We show the 
results in the two more easily trackable limits: one when 



the two hole FSs are completely equivalent, and the other 
when one of the two hole FSs is much weakly coupled 
with electron than the other one and can be neglected. 
In the second limit, 4-pocket model reduces to 3-pocket 
model 28,57 . We show that the system behavior is iden- 
tical in the two limits, and make a conjecture that it 
doesn't evolve between the limits. We continue with our 
earlier assumption of circular electron FSs that nests with 
the hole FS. 

We begin with the limit when one hole FS can be ne- 
glected and 4 pocket model reduces to 3 pocket model. 

A. Effective model with one hole at (0,0). 
The interaction Hamiltonian for the 4-pocket model is 



(i) 



2tt 



Hint — u l ^ c p lS fives' flpts' Cp 3S + U 2 ^ c p 1 sflv 2 s' C Pis'flp 3 s + ^ ( C Pl « C p2 s' flP4.s' flp 3 s + h.C.J 



+ fi <-> fi and «-» uf' 



,(2) 



E_£ t t I "4 f t ft J! f _i_ V "4 ft ft f f 

o C pis C y,s' C P4.s' c p 3 s + / J 2 JlpisJlp 2 s' J^Pi" 1 Jlp 3 s ~r 2_^, 2 •>2p 1 sJ2p2s' J2p 4 s' J2p 3 s 

U S flp 1 sf2p 2 s'f2p4.s'flp 3 s + u 7 flp lS f2p2s'f^PiS'f2p 3 s + — \flp is flp 2 s' fapis' fipss + h.C^j 



(15) 



This is a straightforward generalization of the 2-pocket 
case (see also Ref. I57T ) . The notations are the same as for 
the 2-pocket model, but now /i and fi refer to fermions 
from the two different electron bands. The new terms 
uq,ui, and ug are different inter-pocket interactions be- 
tween /—fermions. Because we now have two different 
sets of electron state, it is convenient to work in the un- 
folded BZ, and in Fig [7] we show the interactions that 
contribute to the pairing vertex. There are, however, 
subtle effects related to the actual, As-induced differ- 
ences between folded and unfolded zones, and we will 
discuss them below. 

The two electron bands are related by symmetry k x o 
k y (i.e., £f 1 (k a; , k y ) — e f 2 (k y ,k x )), and it is natural to set 

v!^ = u$ = Ui (i runs between 1 and 4). We verified 
that no new terms are generated under RG flow, how- 
ever the interactions between electron pockets must be 
included as they anyway are generated by RG. 

The angular dependence of the vertices is incorpo- 
rated in the same manner as described in Sec. HH by 
including a cos 2<p terms into the vertices which involve 
fermions near electron pockets. To simplify calculations, 
we first neglect the angular dependence of the intra- 
pocket electron-electron interaction u 4 . Later we show 
that including angular dependence of U4 will not change 
the results qualitatively. 



1. The Vertices 

The computational procedure is the same as before. 
We introduce infinitesimally small SC and SDW vertices, 
dress them up by the interactions, and express the renor- 
malized vertices in terms of the running couplings. 

The diagrammatic expressions for the renormalized 
vertices are presented in Fig. [S] For the SDW ver- 
tex we have A SDW = Af DW + A% DW cos24>, <j> be- 
ing the angle along the electron FS. For SC vertex on 
the hole FS A sc = Aft and on the two electron FSs 
A sc = A e ± A e cos2(f), as required by symmetry for an 
s— wave gap. If |A e | > |A e |, the SC gap has nodes along 
the FS. The nodes are 'accidental' in the sense that they 
are not protected by any symmetry. 

For the SC vertex we then obtain the coupled set of 
equations 

1 - u 5 L -2u 3 L -2cm 3 L \ / A° h \ I A h \ 
~u 3 L l-u 4 L ) = A e [16) 

-2au 3 L 1 J \ A° J \A e J 

where M4 = 114 + itg. For the SDW vertex we obtain 
(assuming that Ai^ are real) 

/ 1 + («i + u 3 )L f ( Ul + u 3 )L \(A° 1 \_(A 1 X . 
\ a(u 1+ u 3 )L 1 J { A° J - { A 2 } U > 
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FIG. 6: Running vertices in SDW and SC s± channel for the 2-pocket case at perfect nesting as functions of L. Three 
qualitatively different scenario are possible: (a) vertices diverge before Le f is reached. Both vertices flow to the same fixed 
point, their ratio tends to one, and the fixed point has an extended symmetry, (b) Le f is reached before the fixed point, but 
in the range of L where the SC vertex is already attractive (positive). The vertices flow independent of each other beyond 
Le f (below Ef), the SDW vertex diverges first. On doping SDW yields to SC (see Fig. 0). (c) The value Le f is reached when 
SC vertex is still repulsive (negative). In this case SC instability does not occur even after SDW instability is eliminated by 
doping. 




FIG. 7: Hole (center) and electron (edges) FSs in the unfolded 
BZ for 4-pocket model. The arrows with symbols indicate 
various intra-pocket and inter-pocket interactions which con- 
tribute to the SC vertex. There exist other density-density 
and exchange inter-pocket interactions (ui, U2, U6, and u-j 
terms, not shown). 




FIG. 8: The diagrams for the SDW and SC vertices for the 
4-pocket model (panels (a) and (b)). The notations are the 
same as in Fig(4] Single and double dashed lines describe 
fcrmions from the two electron bands. 



s± 



-{v,i + U5) + \J (u4 - ws) 2 + 8(u 3 ) 5 



We included angular dependence of both u 3 and m 
terms and for simplicity set a to be the same for both. 
Observe that a— dependent terms appear in SDW vertex 
with extra 1/2 compared to SC vertex. This is because 
internal and external part of the SDW vertex each con- 
tains one fermion from the electron FS, while in the SC 
vertex there are two such fermions, either in the internal 
or in the external part. One can easily check that in the 
SC vertex all a— terms then appear with extra factor of 
2 compared to SDW vertex. 



The effective vertices are found by diagonalizing these 
matrices and casting the results into the forms Aj = 
A°(l + TiL). For a = 0, the formulas simplify and we 
have 



V s c = T s++ 



-(U4 + U 5 ) - y/(u 4 - U 5 ) 2 + 8(u 3 ) 2 



r 



SDW 



= r? 



SDW 



= U 1 +U 3 



(18) 



The solutions corresponding to T% and T2 are A^/A e < 
and Ah/A e > 0, accordingly, hence the notations T s± 
and r s++ . The vertex r s++ is repulsive for all couplings, 
while r s± is repulsive for 2u\ < U4U5 and is attractive 
for 2u§ > U4U5. The SDW vertex is attractive. We re- 
call that U4 and u§ are Coulomb interactions at small 
momentum transfer, while u 3 and u% are Coulomb inter- 
actions at large momentum transfer and are supposed to 
be smaller than it 5 , u 4 . At the bare level then there is no 
attractive component of T sc for a = 0. 

At a finite a we obtain two SDW vertices rf DW = 
(tti + u 3 )(l + Vl + 2a 2 )/2 and T^ DW = + 
2a 2 — l)/2. The second one is repulsive, which 



« 3 )(vT 
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the first one is attractive and only increases with a. For 
the SC vertex, the three rf are obtained by diagonal- 
izing Eq. (|16j) what requires solving the cubic equation. 
The analytical expressions for Tf are long and we re- 
frain from presenting them (we will show all three Tf c in 
the figures). It is essential, however, that for any a^O, 
one of three Tf c is attractive even when is larger 
than 2u§ and at a = SC vertices are repulsive. At small 
a, we have for such induced solution 



(a) 



SC 



4U3U5 



U4M5 



2«| 



1*41*5 > 2u 3 



(19) 



Other two solutions rf 3 are negative, i.e. repulsive. 

For U4U5 < 2itf, the attractive solution exists already 
at a = 0, and is only weakly affected by a. For 1*41*5 < 
2it§ we have: 



r? c = r 



s± 



1 



'2u2(-Aw+ V(Au) 2 + 8u 2 )' 



8u§ 



(20) 



\ (r s± )V(Au) 2 

where Am = U4 — 1*5 and r s± is given by Eq. [181 The 
her two rf 3 are again negative. 
We now proceed with the RG flow. 



2. RG flow between A and Ef 

Like in 2-pocket case, at A > E > Ep, renormal- 
izations in p-h and p-p channels are logarithmical and 
independent of the location of external momenta. The 
derivation of parquet RG equations is straightforward 
but requires more efforts as there are new terms in the 
Hamiltonian. For illustration, we show in Fig [S] the dia- 
grams contributing to the renormalization of the vertices 
1*4 and 1x3. The diagrams for the renormalization of other 
vertices are similar. 

Collecting the diagrams for the renormalization of all 
couplings, we find that the terms uq ± and 1*4 — ug, 
are decoupled from the rest of the terms and are renor- 



malized as iij 



Because all these uj are the 



differences between Coulomb interactions at small and 
large momentum transfers, their bare values are positive 
in which case the these interactions flow to zero under 
RG and are therefore irrelevant. 

The other five vertices are all coupled and flow accord- 
ing to 



its 

7/4 = 

iii = 

U2 = 

U3 = 



•2«i] 



= - [&l + 2uj] 

[uj + 4] 

[2-U1-U2 — 2u|l 
[Auiu 3 — 2u 2 u 3 



U5-U3 - U4U3J 



(21) 



This set of equations can be easily solved numerically. 
Fig ITUl shows the plot of U5/1/1, U^/ui, 1*2/^1, and 1*3/1*1 
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FIG. 9: Second order diagrams for the renormalizations of U4 
and 113 vertices (panels a and b, respectively). The combina- 
torial factors are not shown but must indeed be included. 



with L. For simplicity, we set bare values of U4 and 115 
to be equal - the two then remain equal in the process of 
RG flow. The values of the ratios at the fixed point are 
indicated by the dots. These can be easily found analyti- 
cally by requesting that all 5 equations in (|21|) be identi- 
cal. Imposing this condition we obtain 1*5/1*1 = u^/ux = 
-y/6 « -2.45, u 3 /ui = (3 + 2V6) 1/2 « 2.81, u 2 /u x = 0. 
Like in 2-pocket model, intra-pocket repulsions 114 and 
its decrease under RG, change sign at some value of 
-L, and become negative at larger L. This sign change 
(overscreening) goes beyond a conventional McMillan- 
Tolmachev screening of the Coulomb interaction, and is 
the result of the "push" from 1*3, which in turn increases 
under RG due to the "push" from u\ which contributes 
to the SDW vertex. So, eventually, overscreening is the 
result of the "attraction" between SDW and SC fluctu- 
ations. The difference U4U5 — 2u 3 also changes sign at 
some L and becomes negative at larger L. 

We now substitute the running couplings into the ex- 
pressions for SC and SDW vertices and check how they 
flow. The results are presented in Fig. [TT] for a = 0.4. 
For the SDW vertices, the positive one increases with L, 
like in the 2-pocket model, while the negative one be- 
comes more negative, i.e., even less relevant. For the SC 
vertices, Tf c is positive for all L, the other two rf 3 are 
negative and hence irrelevant. The positive Tf c interpo 



lates between Eq. (fl"§)) at small L, and Eq. (f2T)]) at larger 
L. We emphasize that for all values of L this is the same 
solution, i.e., there is no level crossing (see Fig. [TTI) . 

In Fig. [T2k we compare the behavior of T SDW and T sc 
as functions of L assuming that the fixed point is reached 
at energies above Ep. At small L, we have the same sit- 
uation as in 2-pocket model: SDW vertex is larger than 
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FIG. 10: The solution of Eq. (gTJl - the RG flow of u 5 , u 4 , 
U2, and U3 (all relative to Ui) with L. For simplicity, we set 
bare values of Ua and us equal. Note how, as fixed point is 
approached, the renormalized Coulomb repulsion at small mo- 
menta (us and in terms) is suppressed and eventually changes 
sign, while the pair hopping term us is strengthened. 



vertex. For the value of a which we used in this figure 
(a = 0.4) the transition from nodal to nodeless gap oc- 
curs at L smaller than the one at which SC vertex crosses 
the SDW vertex, i.e., when SC becomes the leading in- 
stability, the superconducting gap is already nodeless. 
But for other values of a we can get either nodeless or 
nodal SC in this regime. In Fig. [T3]we plot the "phase 
diagram" coming out of parquet RG for different a [the 
bare values of u, are the same for all figures in this sec- 
tion]. In white region, SDW vertex is the largest and 
SC vertex is subleading, implying that superconductiv- 
ity can be revealed only after SDW order is suppressed 
by doping. In the shaded region the SC vertex is the 
largest. We see that superconducting gap in this region 
can actually be either nodal or nodeless, depending on 
the value of a. At a — 0, the ratio of A e and at the 
fixed point is A/j = — V2A e . 



3. RG flow below the scale of Ef 



SC vertex. However, the rate of increase of the SC ver- 
tex exceeds that of the SDW vertex, and at some L be- 
fore the fixed point is reached, T sc crosses over Y SDW 
implying that superconductivity becomes the leading in- 
stability even at perfect nesting. Such crossing has been 
reported in fRG calculations^ for the same model. We 
view the agreement as a good indication that numeri- 
cal fRG and analytical parquet RG approaches describe 
the same physics. In our analytical RG, the reason for 
the crossing is combinatoric: compared to 2-pocket case 
(where SDW and SC vertices flow to the same value un- 
der RG) , the presence of the second electron FS adds the 
factor of 2 to the renormalization of the SC vertex as a 
pair of fermions from the hole FS can hope to each of 
the two electron FSs. However, there is no such factor 
of 2 in the renormalization of the SDW vertex due to 
momentum conservation. 

We emphasize that the crossing of rf DW and Ff is 
not related to the angular dependence of the interaction. 
Even when a — 0, SC vertex exceeds SDW vertex near 
the fixed point of parquet RG. At the fixed point, the 
ratio of the two is rf c /Tf DW 1.69. 

The SC order parameter by itself has an interesting 
character. We recall that we approximate the gap along 
the hole FS by a constant A^ and approximate the gap 
along the two electron FSs by A e ± A e cos 2(f). At small 
L, the attractive T sc exists only because of a non-zero 
a, and A e is larger that A e (Ref. |28I ). hence the gap 
along the two electron FSs has "accidental" nodes. As 
L increases, the SC vertex T evolves, according to Fig. 
[TTl and eventually gets close to the would-be solution for 
a = 0. For the latter, A e = 0, and the gap obviously 
has no nodes. The crossover from one limit to the other 
is displayed in Fig. [T3J where we show the flow of the 
gaps Ah,A e , and A e , corresponding to the leading SC 



We now consider the situation below Ep . As in the 
2-pocket model, we have to introduce three different 113 



couplings {u^\ 



lb) 



3 u; , and (U3 ) of which u 3 u> is the part 



of SDW vertex, and tig is the part of the SC vertex (the 

corresponding replace u 3 in Eqs. (jTI)]) and (fT7|). The 
flow of the couplings is now governed by 



«5 = - [u§ + 2(4 C) ) 2 

u 4 = - u\ + 2(4 c) ) 2 

u x = + [ul + (uf ) ) 2 

ii 2 = + [2uiu 2 — 2ul] 

U3 = 2uiu^ — 2u2U^ 

■ (b) o ( b ) 

U 3 — ZM1M3 

• ( c ) r 1 — 1 (c) 



(22) 



Note the M4 and 1*5 have identical equations and hance 
treated identically. One can then straightforwardly ver- 
ify using (|2"2"j) that SC and SDW vertices decouple, as 
they should, and each satisfies T; = rf. Hence, as be- 
fore, whichever vertex is larger at Ep gives rise to the 
first instability as T decreases. If SC vertex prevails, the 
system becomes SC at perfect nesting and remains a SC 
at finite dopings (Fig. IT2"b). If SDW vertex prevails, 
the system becomes an SDW antiferromagnet at perfect 
nesting and then eventually becomes a SC upon doping 
(Fig. [T2"b). In distinction to the 2-pocket case, we don't 
need to worry about the sign of the SC vertex once SDW 
instability is reduced by doping because one of T sc is al- 
ways attractive (see Fig. [TT| . We emphasize again that 
this attractive T sc leads to either nodeless s± gap, or 
to s± gap with nodes along the electron FSs, depending 
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FIG. 11: The flow of SC and SDW vertices under RG in the effective 3-pocket model for a — 0.4. Panel (a) 3 SC vertices. 
One solution is attractive for all L (and corresponds to s± pairing), the other two are repulsive. One of the repulsive solutions 
is of s± character, another of s + + character. At small L, the positive solution is the one induced by a, at large L it almost 
coincides with the solution which becomes positive for these L already at a = 0. The circle marks the area where positive and 
negative solutions come close to each other. The splitting between the two increases with a. (b) The two SDW vertices. One 
vertex is always attractive (positive) and the other is repulsive (negative) . 
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FIG. 12: The flow of vertices for different values of Le f for a = 0.4. (a) RG flow reaches fixed point before Le f - SDW vertex 
is larger at small L, but SC vertex 'crosses' the SDW vertex at some distance from the fixed point and becomes the strongest 
vertex at the fixed point. As a result, the system develops a SC order, (b) Le f is reached after the 'crossing' but before 
reaching fixed point. The system still develops a SC order, and SDW order does not emerge, (c) Le f is reached before the 
'crossing'. In this case the SDW vertex still develops at small dopings, and SC order emerges at larger dopings, when SDW 
order gets suppressed. 



on a and on the interplay between Ep and the scale at 
which parquet RG flow reaches the fixed point. 



tion U4 = 114 + ug, 
We have 



so we need to consider only this term. 



Effect of the angular dependence of electron-electron 
interaction 



For completeness, we present the results for the evo- 
lution of the SC gap structure under RG flow for the 
case when we preserve the angular dependence in the 
electron-electron interactions - 1*4 and us terms. These 
interactions only contribute to the pairing channel, so it 
will be sufficient co consider U4 and u% interactions be- 
tween fermions with momenta k, —k;p, —p. The generic 
structure of the angular dependence of such interactions 
is given by Eq. We found earlier that M4 and ug 

terms contribute to the s— wave pairing in the combina- 



Ui{k,p) = ua (l + 2cv' (cos20/c ± cos2c/) p ) ± 4a cos2</>/c cos 

where plus sign is for intra-pocket interaction and minus 
sign is for inter-pocket interaction. 

There are two new effects associated with the angu- 
lar dependence of u^{k,p). First, when U4U5 > 2u§, the 
pairing vertex T sc not necessarily has an attractive com- 
ponent. It was always the case for angle- independent U4. 
Now the existence of the attractive interaction is subject 
to condition K > 0, where 



K = 2u 4 u 5 ((a') 2 - a") + ul{a 2 + 2a" - 3aa'). (24) 
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FIG. 13: The evolution of the three components of the super- 
conducting gap with L for a = 0.4. is the gap on the hole 
FS, and A e ± A e cos 2<f> are the gaps along electron FSs. A 
circle marks the point where A e and A e cross, and the gap 
along each of electron FSs changes from nodal to nodeless. 
Note that the solution is always of s± character; meaning A^ 
and Al axe of opposite signs. 
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FIG. 14: The RG flow of the 4-pocket model in variables a 
and L over the range of energies above Ef- The white zone 
is where SDW vertex is the largest, and the shaded zone is 
where SC s± vertex is the largest. The thick solid line in the 
shaded zone marks the transition form nodal to nodeless s± 
gap in the region where SC vertex wins. The dashed line is 
the continuation of this transition line into the region where 
SDW vertex takes over. 



If we set the angular dependence of 1*3 and W4 terms to 
be equal, i.e., set a' = a, a" = 0, this condition reduces 
to > 1*3 , which is well satisfied. However, for a 

generic a' and a", Eq. (j2~4")l is not necessarily satisfied, 
and UK < 0, attractive Tf c appears only above some 
RG scale L, like in 2-pocket model. 

Second, the gap structure may change in some range 
of L. To demonstrate this, make angular dependence of 
U3 and U4 equal, i.e., set a' = a, a" = 0. The set of 
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FIG. 15: The phase diagram of the 4-pocket model at perfect 
nesting in variables a and Le f = logA/_B_F- In red is the 
region where SDW order develops, and in green is the region 
of the SC order. Dark and light green regions correspond to 
nodal (dark) and nodeless (light) SC gap along the electron 
FSs. The two sub-figures show the behavior at finite doping. 
Superconducting state which is brought out upon doping in 
the left sub-figure is either nodeless or nodal depending on 
the location with respect a dashed line inside the red (SDW) 
region. The transition between SDW and SC states can be ei- 
ther first order or involve intermediate co-existence phas o 58 ' 59 . 



equations for the SC vertices then becomes 



1 — u^L —2u^L — 2au^L 
—113L 1 — U4L — CXU4L 
—2au^,L —2auiL 1 




A. 



25) 



As before, we need to diagonalize this set, cast the result 
in the form A, = A°(l-rTjZ/) and consider the largest IV 
The evolution with L of Aft, A e , and A e for such is 
shown in Fig 1161 and the phase diagram is shown in Fig. 
[T71 We see that, over some range when A e is the largest 
and the gap has nodes along the electron FSs, the gap 
actually has "nodal s + +" character in the sense that 
Aft and A e are of the same sign, although the dominant 
term is still the oscillating component A e . Note, however, 
that the character of the gap changes back to s± before 
it becomes nodeless. 

This appearance of the nodal s + + like gap might seem 
unusual, but it should be kept in mind that this gap is 
present in the parameter range where without the angu- 
lar dependence there wouldn't have been a solution. The 
firm requirement then is that in the solution induced by a 
the oscillating A e component is the largest, because this 
is the way to minimize the effect of intra-pocket Coulomb 
repulsion. The relative sign between the subleading Aft 
and A e terms is not uniquely determined by this require- 
ment and be either minus or plus, depending on the in- 
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FIG. 16: Behavior of the different A's when the angular de- 
pendence of the electron intra-pocket coupling U4 is included. 
As before, we set a = 0.4. The only difference compared to 
Fig. [I3]is the appearance of the region, at small L, where SC 
order parameter has nodal s + + character meaning that Ah 
and A e are of the same sign. In this range of L the SC vertex 
is, however, smaller than the SDW vertex. The character of 
the SC gap changes to nodal s± and then to no-nodal s± 
before SC vertex takes over SDW vertex. 
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FIG. 17: Same as in Fig. 1141 but with the angular depen- 
dence of the electron intra-pocket coupling U4 included. The 
only new feature is the existence of a range where SC vertex 
(secondary in this range to SDW vertex) has nodal s++ char- 
acter meaning that the gaps along electron FSs have nodes, 
but the average value of the gap along the electron FS is of 
the same sign as the gap along the hole FS. 



B. The 4-pocket model in the limit when two hole 
FSs are identical 



We now consider the opposite limit where the two hole 
FSs centered at (0,0) are equivalent. We show that the 
system behavior in this second limit is the same as in the 
first. The equivalence of the two limits hints that the 
system behavior in the intermediate case is very likely 
the same as in the two limits. 

The computations in the case of two identical hole 
FSs proceed in the same way as before, but there are 
more vertices. The new terms are u\, U2, and u 3 interac- 
tions with hole fermions from the second hole bands, the 
analogs of these three interactions for fermions from the 
two hole bands, U5 interaction for the second hole band, 
and the interactions of the kind c\c\ciC2- Note that 
there are no /1/1/1/2 terms for fermions from the two 
electron bands because they would violate momentum 
conservation. 

The full set of RG equations is rather cumbersome, 
but we verified that (i) RG flow indeed preserved the 
invariance between the two hole bands, and (ii) all intra- 
pocket and inter-pocket interactions involving fermions 
from the hole bands flow to the same value U5. The 
analysis based on 5-orbital Hubbard model also yields 
near-equivalence of all it, involving fermions from hole 
pockets^. 

To simplify the presentation we set all interactions in- 
volving fermions near hole FSs to be equal to u 5 from 
the start and also neglect the angular dependence of the 
interactions. We also set M4 = u$ because u 4 — ug > 
again flows to zero under RG (see paragraph before Eq. 

EE}. 



1. The Vertices 

The equations for the SC and SDW vertices are ob- 
tained in the same way as before (see Fig. ©), but now 
A/i is composed from fermions with k and — k belonging 
to either of the two hole pockets. This leads to the equa- 
tions for the SDW vertex Ai and SC vertices A^ and A e 
in the form 



Ai = A?(l+Ui+U 3 ) 



(26) 



and 



terplay between electron-hole and electron-electron inter- 
actions. 

These two potential changes introduced by the angu- 
lar dependence of U4, however, affect only the behavior 
at small L. At large L (i.e., at low energies), the system 
behavior remains unchanged: SDW vertex is the largest 
at small/intermediate L, but SC vertex still crosses over 
SDW vertex at some L, and beyond this scale SC insta- 
bility comes first. 



1 - u 5 L ~u 3 L \{A° h 
-2u 3 L l-u 4 L \ A° 



A h 
A R 



(27) 



where u\ = 2u\, u 3 = 2u 3l 114 = 2u 4 and U5 = A1L5. 
Casting the results into Aj = A°(l + I\L), we obtain 



T SDW = u 1+ u 3 



-(u 4 + «s) + ~ U5) 2 + 8(w 3 ) 2 



(28) 
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2. Parquet RG equations 

The RG equations are obtained in the same way as 
before and are 



115 = 

U3 



L~,2 



+ [4{tiu 3 - u 5 {t 3 - W4U3] 



(29) 



We drop ui for simplicity as it eventually becomes 
smaller than other m 

Comparing these equations and the equations for the 
vertices with those for one hole FS (Eqs. [T5] and |2"T|) 
we see that they are identical up to renormalizations 
Ui — > iii. Accordingly, the flow of the couplings and the 
vertices is the same as in the limit when only one hole 
FS is present. In both limits, SC vertex is secondary to 
SDW vertex at large energies, but has larger slope and 
crosses over SDW vertex at some energy, before the sys- 
tem reaches a fixed point. At smaller energies, SC vertex 
is larger, i.e. if parquet RG flow extends beyond the 
scale where the two vertices cross, the system first devel- 
ops a SC order even at perfect nesting. This SC order 
can be either with or without nodes in the gaps along 
the two electron FSs (see Fig. [T5]) . The only difference 
to the effective 3-pocket model is that now at the fixed 
point we have A e = — V2Ah for a = as opposed to 
Aft = — V^Ae for the earlier case. 

As we said, the equivalence of the system behavior in 
the two limits strongly suggests that the same behavior 
holds also in the intermediate case. 



SDW dominated region (red), SC emerges after doping 
reduces SDW order. In the other part where SC domi- 
nates, SC order prevails already at perfect nesting. The 
s± SC gap can be nodeless or have nodes along electron 
FSs depending on how strong is the angular dependence 
of the interaction between electrons and holes. 

A final remark: In the analysis above we considered 
only the interactions which obey momentum conserva- 
tion in the unfolded BZ. These are direct interactions 
between fermionic states obtained by the hybridization 
of 5 Fe orbitals. There exists, however, additional inter- 
actions which involve pnictide orbitals as intermediate 
states. These additional interactions obey momentum 
conservation in the folded BZ, but they do not always 
obey momentum conservation in the unfolded, Fe— only 
BZ. An example of such process is shown in Fig. [T5J 
two fermions from the hole band near k = scatter into 
two fermions at two different electron pockets. In the 
unfolded zone, this process doesn't conserve momentum, 
and we didn't include it into our consideration. In the 
folded zone, both electron FSs are at (tt, tt), and this pro- 
cess is an umklapp process. The difference is indeed due 
to the fact that in reality such process involves interme- 
diate states on As. 

Neither our RG procedure nor fRG calculations include 
such terms. How important are they is not known. On 
general grounds, such interactions tend to enhance the 
SDW vertex and might potentially alter the picture that 
we presented. They also may alter the ordering momen- 
tum of the SDW state. This remains an open issue. 



C. Summary of the results for the 4-pocket model 

Collecting all the points we have discussed - we have 
shown that under suitable extent of renormalization of 
the Coulomb repulsion and pair-hopping couplings one 
can have SDW, nodal s±, and nodeless s± state even 
at perfect nesting. The angular dependence of the in- 
teraction between holes and electrons tends to drive the 
system towards a nodal SC phase. The SC order devel- 
ops if the fixed point is reached within parquet RG, but if 
the scale Ep is reached before that, the system develops 
either SDW or SC order (either nodeless or nodal), de- 
pending on at what L the flow crosses over from parquet 
to ladder RG. That the SC s± order can emerge even at 
perfect nesting is specific feature of the 4-pocket model. 
This feature was not present in the 2-pocket model, where 
the fixed point had an 0(6) symmetry. This symmetry is 
clearly broken in the 4-pocket model, even when a = 0. 
The 'crossing' of the SDW and SC vertices can be unam- 
biguously attributed to the presence of the other electron 
pocket because its presence helps SC but not SDW. 

Fig. [15] summarizes the implication of our results to- 
wards the actual phase diagrams of Fe-pnictides. In the 



FIG. 18: The scattering which takes two fermions from the 
region near k = and scatters them to fermions with mo- 
menta near (7r,0) and (0, tt). This process is not allowed in 
the unfolded BZ because of momentum non-conservation, but 
it is allowed as an umklapp process in the folded BZ, which 
knows about As. 



V. 5-POCKET MODEL 

We now extend the analysis from previous two sections 
to a 5-pocket model in which we include into considera- 
tion the additional hole pocket appearing at (tt, tt) point 
in the unfolded BZ. We show below that the behavior of 
5— pocket model is similar to that for 2-pocket model in 
the sense that SDW vertex exceeds SC vertex along the 
whole RG trajectory, and SDW and SC vertices tend to 
the same value if the fixed point is reached within parquet 
RG. 
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As in the previous section, we restrict our consideration 
to the two limits, one when the two hole FSs centered 
at (0, 0) are identical, and the other when one of these 
two hole FSs is relatively weakly coupled to electronic 
states and can be neglected. In the latter case, 5-pocket 
model reduces to an effective 4-pocket model consisting 
of one hole FS at (0,0), one hole FS at (tt,tt), and the 
two electron FSs at (0,tt) and (it, 0). We show that the 
system behavior is again identical in the two limits. 




FIG. 19: The FSs and interactions in the 5-pocket model. 
Dashed lines mark the interactions already present in 4-pocket 
model, solid lines mark the new pairing interactions specific 
to 5-pocket model. As before, we only present interactions 
which contribute to the pairing vertices. There are other 
density-density and exchange interactions between electrons 
belonging to different pockets. 



A. Effective model with one hole pocket at (0, 0) 



The FS geometry and interactions contributing to the 
SC vertex for the effective 4-pocket model with hole pock- 
ets at (0, 0) and (ir, n) are presented in Fig. [19] 

The Hamiltonian now contains three new terms itg, 
Mio, and itn, which are density-density, exchange, and 
pair-hopping interaction between fermions belonging to 
two different hole pockets. In addition, we have three 
new vertices shown in Fig. 1201 These include fermions 
from two different hole and two different electron FSs. 
We call them Wi vertices (i runs from 1 to 3). 




FIG. 20: The new interaction vertices for the 5-pocket model. 
Single and double solid lines denote fermions from the two 
hole pockets, single and double dashed lines denote fermions 
from the two electron pockets. 



The Hamiltonian now has the form 



11 int 

Zn 



(i) 

= U l ^ C lpi sflvos'hpAs'Clpzs + u\ ^ c| pi s /i D9S 'Clp 4 s / flp 3 s + f c l Pl s C \v^s' f^iS' Ilp 3 s + h.i 

+ ft O fi (with c\ unchanged) and u\ ■<-> u[ 2 ^ + c\ «-» C2 (with ft unchanged) and u\ O u\ 
+ cj f>C2 (with fi unchanged) and uf^ O uf^ 

(1) (2) 

+ ^ t o f lpi S f lp 2 S' /lP4 s'flp3S + ^ t Q f2pxsf2p 2 s' f 2 P&s' f 2 P3S 



2 



,(2) 



E 5 t t \ ^ 

~^~ c lp 1 s c lp 2 s' c ^P4.s' c lp 3 s + 2 , ~~^~ c 2p l s c 2p 2 s' c 2piS l c 2p 3 s 

+ Ue flpisf2p 2 s'f2p4.s'flp3S + u l fl Pl sf2p2s'fipts'f2p 3 s + (^fl Pl sflp 2 s'f 2 P4.s'f2p 3 s 

+ c <-» / and (u 6 , u 7 , u s ) -H> (u 9 , u w , u n ) 

+ J2 Wl C lpisflp 2S 'hp 4S 'C 2 p 3 s + J2 W2 C lpisflp 2S ' C 2 Pi s'f2p 3 s C lpi s4p 2S 'hp iS 'f 2 P3 s 2\30) 

I 



The analysis of the 5-pocket model parallels that of the 
4-pocket model, so we will be rather brief and present 
only the results. We verified that RG equations for 
M4,U5,U8 and un are identical, and these four couplings 
tend to the same value at the fixed point. To make 
the presentation compact, we set U4 = U5 — ug, = Utt 



from the start and call all of them U4,. Similarly we set 
uq = ur — uq — Uto calling it u§. It is convenient to 
introduce = Ui + Wi (i = 1 — 3), and Hi — Ui~ Wi. We 
will use these variables below. 
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1. The Vertices 

We first consider the case when the interactions are 
angle-independent (a — 0) and then discuss system be- 
havior at a nonzero a. 

The SC and SDW vertices are obtained in the same 
way as before, but there are additional terms for the SDW 
term due to to, vertices (see Fig. l2Tj) . Combining this 
with the equations for the SC vertices at (0, 0) and (tt, tt) 
(A/,.1 and A^2 respectively) we obtain 



1 — 2u^L — 113L 
—IU3L 1 — U4L 



-2u 3 L 



-u^L 1 — U4L I 





are attractive, the other two are saddle points. Which 
fixed point the system will flow to depends on the ini- 
tial conditions. In our case all interactions are positive 
(repulsive), i.e. at the bare level Ui are all positive and 
Ui > Ui. For these initial conditions, we verified that 
the system is outside the base of attraction of the second 
fixed point as it can be reached only if bare Wi are neg- 
ative (at this fixed point in vanishes compared to iii i.e., 
Wi/ui tends to —1). 

At the first attractive fixed point Ui vanishes compared 
to Ui i.e., Wi/ui tends to 1. This is consistent with our 
initial conditions. Near this fixed point, U3 can be ne- 
glected compared to U3 in the equations for 114 and uq, 
and the first five RG equations form a closed set: 



and 



Ai = A?(l + («i + u 3 )L) 



(32) 



(we recall that we set u 4 = u§). Casting the equations 
for the SC vertices in the form A, = A°(l + Ti) and 
neglecting repulsive vertex for s ++ SC, we obtain 



T SDW = fii . 

r s± = -2m 4 



"3 
- U 3 



"3 



(33) 



Note that SDW and SC vertices contain different terms 
involving 113. The gap structure for the SC vertex r s± is 
A M = A,, 2 = -A e . 



2. RG flow between A and Ef 

The set of parquet RG equations is obtained in the 
same way as before. Collecting the equations for the 
other variables we obtain 



Ul 
«2 
U3 

2u 6 

2{t4 
U3 
Ul 
f/ 2 



2v,iU2 — 2u\ 

-(2u 6 ) 2 - {u 3 - 
-(2u 4 ) 2 - (fig - 

4uiM3 — 2U2U3 



2m 3 (u 4 + "6) - 2w 3 (w 4 - u e ) 



U3? 

~U3? 



2{t3(u4 



r,2 



-2 



z z =2 

2u\U2 — 2u 2 



■ ue) - 2w 3 (w 4 - ue) 
(34) 



This set of equations almost decouples between the 
subsets for u\ and Ui , the only places where the two sub- 
sets mix are the equations for the flow of u<± and uq whose 
r.h.s. contains both U3 and U3. Re- writing this set as 
equations for the ratios of the couplings, we found four 
fixed points. One corresponds to Ui vanishing compared 
to Ui, another to Ui vanishing compared to Ui, the third 
corresponds to 113 = u 3 , u\ — u\, and the fourth corre- 
sponds to 113 — — M3, u% = ui. The first two fixed points 



2u\ 



U\ — u<i ~r (^3 
U2 = 2u\U2 

ii 3 = 4uiu 3 - 2u 2 U3 - 2it 3 (-u 4 + ue) 
2ii4 = — (2U4) 2 — U3 
2iia = — (2uq) 2 — U3 

Within the same approximation 



YSDW 

ps± 



= Ul - 

-2ui 



U3 
- u 3 



(35) 



(36) 



Comparing these equations with the ones we obtained 
for the 2-pocket model, Eqs. (fTUj) and (fTTj) . we see that 
they are equivalent, up to overall renormalizations of the 
couplings, if we identify 2uq in Eq. (|35[) with U5 in Eq. 
(fTTj) . There is minor difference between the T s± in the 
two cases (2«4 in (|36[) vs U4 + U5 in[T0|). but it vanishes at 
the fixed point. Acordingly, the RG flow of the couplings 
and the vertices is the same as in the 2-pocket model, 
namely SDW vertex remains dominant for all L up to 
a fixed point, and SC vertex changes sign at some L, 
becomes attractive at larger L and becomes equal to the 
SDW vertex at the fixed point if, indeed, the fixed point 
is reached within parquet RG. This similarity with a 2- 
pocket model was first noted by K. Hauled and can be 
understood if we note that in the 5 pocket case SC pairing 
is the same as in 4-pocket model (with extra combinatoric 
factor of 2 compared to 2-pocket case) , but SDW pairing 
is now possible between the two sets of electron pockets 
(see Fig. [2~T|) . this adds the combinatoric factor of 2 also 
to the renormalization of the SDW vertex. 

We now need to understand what is the basis of at- 
traction for this fixed point. For this, we consider the 
two other fixed points for which u 3 = w 3 , u\ = u±, or 
U3 = —U3, ui = ui. We show that both are saddle 
points, and both are unstable when the bare V4 and Wi 
are all positive. 

Consider for example the fixed point at ui = ui and 
U3 = U3. At this fixed point u±/ui = —3, U3/U1 = VTE 
and Uq — 0. Expanding in 5 = Ui — «i and e = u 3 — u 3 , 
we obtain the set of coupled linear differential equations 
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S = 2iii (s + VT5, 
e = 4ui f >/l5<5 + 4e 



(37) 



together with Ui = >uf + m| = 16iif. The solution of 
the set is e = eo ("i) 7 , S = So (ui) J . Substituting and 
solving for the set of two linear equations for eo and So, 
we obtain 71 = 11/8 and 72 = —1/4. For e and 6 which 
correspond to 7 = 71, the fixed point is unstable, for 
7 = 72, it is stable. A simple analysis shows that 71 
is the solution when eo/<5o > 0, while 72 is the solution 
when €q/Sq < 0. In our case, the bare values of e and S 
are 2ws and 2wi, respectively, both are positive. Hence 
this fixed point is unstable, and the RG flow bring the 
system towards the stable fixed point at which u\ and 113 
are both small. The stability analysis of the fixed point 
at u\ = ui, or U3 = — U3 yields the same results, leaving 
the fixed point with Ui << as the only stable fixed 
point. 

We next consider how the results are modified due to 
angular dependence of the vertices. We found two effects. 
First, one of SC vertices Tf c can become attractive al- 
ready at small L in the same way as in the 3-pocket 
model studies in the previous section. Namely, the sys- 
tem adjusts cos2(/> and angle-independent components of 
the gaps along the two electron FSs to minimize the effect 
of intra-pocket Coulomb repulsion. Just as for 3-pocket 
model, rf c is attractive and scales as a 2 if we only in- 
clude angular dependence of the pair-hopping 1*3 and W3 
terms. Second, SC and SDW vertices do not become 
identical at the fixed point if a is nonzero. If we only in- 
clude angular dependence of U3,u>3 and u\ and w% (and 
set them equal), we find that SC vertex becomes larger 
than SDW vertex very near fixed point. However, the 
effect is numerically very weak, even when a ~ 1. In 
Fig. [2U we show the flow of SC and SDW vertices for 
a = 0.3. SC vertex eventually becomes larger, but this 
is truly weak effect. 

The flow of SDW and SC vertices towards almost the 
same value in the 5-pocket model has been found nu- 
merically by Thomale et al within fRG study2£. Our 
analytical RG results for this case again agree with their 
fRG, what, in our opinion, is another confirmation that 
the "topology" of the RG flow is chiefly determined by 
combinatoric effects. 



3. RG flow below the scale of Ef 

The system behavior for the case when the fixed point 
of the functional RG is not reached at E > Ep is quite 
similar to the 2-pocket model with the only difference 
that now SC instability always occurs when SDW order 
is destroyed by doping. Namely, at perfect nesting the 
system develops SDW order. At finite doping, the RG 
flow of the SDW vertex levels off, and SC vertex eventu- 
ally becomes larger. The SC gap has s± structure either 



u 1 



w. 




n >4 



FIG. 21: The diagrammatic equation for the renormalized 
SDW vertex in the 5-pocket model. Comparing with the cor- 
responding Fig. [8^ for the 4-pocket case, there are extra 
diagrams which contribute to the SDW vertex. This leads to 
effectively doubling r SDW . 



without or with nodes along the electron FS, depending 
on the values of a and of log A/Ep. The phase diagram 
is similar to that in Fig. [T5J but only has the "SDW" 
region in that figure. 
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FIG. 22: Parquet RG flow of SC and SDW vertices for the 
5-pocket model at a — 0.3. The SDW vertex remains the 
largest over the whole flow, and the ratio of the SDW and 
SC vertices approaches one at the fixed point of parquet RG. 
This is very similar to the 2-pocket case except that here r sc 
is attractive for all L. 



B. 5-pocket model with two equivalent hole FS at 

(0,0) 

We now consider the opposite limit of 5-pocket model 
when the two hole pockets centered at (0, 0) are equiva- 
lent. Our goal is to verify whether the system behavior 
remains the same as in the limit when we keep only one 
of these two hole pockets. 

The computations in the case of two equivalent hole 
pockets at (0, 0) are quite involved and we only present 
the results for a = 0. Because the pockets at (0, 0) and 
(tt, 7r) are now non-equivalent (in the sense that there 
are two pockets at (0,0) and only one at (n, 7r)), the 
interactions involving these pockets do not need to flow 
to the same value under RG, e.g, Ui and W4 need not to 
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FIG. 24: The two non-equivalent SDW vertices in the 5- 
pocket model with two identical hole FSs at (0, 0). The triple 
solid line stands for the fermion (hole) at (tt, it). These two 
vertices are present also for the case when there is only one 
hole FSs at (0, 0), but in that case we set A° = A£ and verified 
that the equivalence also holds for renormalized Ai,2- 



FIG. 23: Comparison of the SC phases for 4 and 5 pocket 
models. The initial conditions are identical for the common 
parameters. The figure shows that there exists a parameter 
range where the SC gap is nodeless in the 5-pocket model but 
is nodal in the 4-pocket model. 



fly to the same value, and also electron-hole and hole- 
hole interactions involving fermions from near (0,0) and 
(it, tt) need not to be the same. Finally, SC gaps on the 
hole FSs at (0, 0) and (tt, tt) also do not have to be equal. 



1. The Vertices 

Keeping all this in mind and applying the same anal- 
ysis as before we obtain the equations for the SDW and 
SC vertices. For the SDW, we introduce two vertices Ai 
and A2 , shown in Fig. [24] and write the set of 2 x 2 
coupled equations as 



l + 2(ui+u 3 )L (wi+w 3 )L 
2(w 1 +wz)L 1 + (ui+m 3 )L 



A? 
A? 



A 2 



38) 



where the vertices u\ and M3 are shown in Fig. 1251 

For the SC vertex we introduce, as before, A e , A/ji 
A(fc = (0, 0)) and A h2 = A(k = (tt, tt)) and obtain 



1 - 2u 4 L - 
-2u 3 L 1 
-2u 3 L 



Au 3 L 
- Au 5 L 
Au 5 L 





The vertices M5 and M5 are shown in Fig[ 



2. Parquet RG equations 

The total number of RG equations is quite large and 
we refrain from writing all of them. Quite predictably, 
the combinatoric factors associated with the existence of 
the two equivalent pockets at (0, 0) give rise to relations 
Tii = 2ui (i = 1,3,5), U4 — 4it5, and U5 = 4«5. Us- 
ing these relations, we obtain for the relevant couplings 
Mi, toi, U3, w 3 and U4 the set 
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FIG. 25: The interactions involving fermions near the hole 
FS at (7r,7r) (the triple solid line). The dashed line stands 
for a fermion from an electronic pocket. For the case of only 
one hole FSs at (0,0) we set us = us = us, u\ = m, and 
^3 = ^3 from the beginning and verified that these relations 
hold for running couplings. For the case of two identical hole 
pockets at (0,0), the fixed point values of the couplings in- 
volving fermions near (tt, tt) are different from those involving 
fermions near (0,0). 
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(40) 



Introducing now U\ = 2u\ 4 
u 3 = 2u 3 + V2w 3l u 3 = 2u 3 
into ([4"U| we obtain 



\/2wi, mi = 2mi — y/2wi, 
— V2w 3 and substituting 
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This is exactly the same set as Eq. ([34)) that we ob- 
tained in the previous subsection (we skip the equation 
on ui which is irrelevant coupling anyway). 

Under the same conditions ( Ui — 2ui (i — 1,3,5), 
u 4 = 4u 5 , and u 5 = 4u 5 ), the relevant SDW and SC 
vertices become 

r SDW = u 1 + u 3 

T s± = -2u 4 + M3 + % (42) 

These again are the same equations as Eqs (|33|) for 
the case of only one hole FS at (0,0). The only differ- 
ence with the other limit is that now the solutions cor- 
responding to T SDW and T s± from (g^ are A 2 = V%Ai 
and A M = 2A hl = -A e . 

We see therefore that the system behavior in the two 
limits is identical. Like in the 4-pocket case, this equiva- 
lence strongly suggests that the same behavior holds also 
in the intermediate, most generic 5-pocket model, when 
the two hole pockets at (0,0) are both present but are 
not identical. 



C. Summary of the results for 5-pocket model 

We see that the system behavior in a 5-pocket model 
is "intermediate" between 2-pocket and 4-pocket mod- 
els. On one hand, like in a 4-pocket model, the largest 
SC vertex can be positive already at the smallest L, even 
when intra-pocket Coulomb repulsion is the dominant in- 
teraction. If this is the case, then there is no critical L 
before which SC vertex is repulsive, and the system al- 
ways becomes a SC when the competing SDW instability 
is reduced. The SC gap is either nodeless or with nodes 
on electron FSs, depending on a, much like in the un- 
shaded region in Fig. [14] On the other hand, like in 
a 2-pocket model, SDW vertex remains larger than SC 
vertex for all L before the fixed point is reached, and the 
two vertices flow to the same value at the fixed point of 
parquet RG (Fig. [22]) . This last statement is exact when 
a — and remains numerically quite accurate even when 
a ^ although strictly speaking, at a finite a, SC ver- 
tex eventually becomes larger than SDW vertex in the 
immediate vicinity of the fixed point. 

D. Comparison of 4-pocket and 5-pocket models 

It is instructive to compare the structures of the SC 
gaps in 5-pocket and 4-pocket models for the same val- 
ues of input parameters (and using the same relations as 
above for extra parameters of a 5-pocket model). This 
comparison is shown in Fig. [331 We see that there is 
quite wide parameter range where in the 4-pocket model 
the gap has nodes while in the 5 pocket model it is still 
nodeless. Each point in the phase diagram in Fig. [521 
corresponds to some values of the couplings, hence this 
result implies that for a certain range of input parameters 



4-pocket model yields a gap with nodes while 5-pocket 
model yields the gap without nodes. This agrees with 
the number of RPA studie o 20 ' 29 ' 30 which found nodal gap 
for 5-pocket model and no-nodal gap for 4 pocket model. 
At the same time, our results show that in both models 
there are regions of parameters in which the SC gap is 
either no-nodal or has nodes. 



VI. CONCLUSION 

We have done calculations addressing on the same foot- 
ing the issues of the interplay between intra and inter- 
pocket Coulomb repulsion in the Fe-based superconduc- 
tors, the competition between SC and SDW orders, and 
the angular dependence of the SC gap. We considered 
2-, 4-, and 5-pocket models for the pnictides and for each 
model considered the flow of the couplings and of SDW 
and SC vertices within analytical parquet RG scheme. 
We found that in all models, fluctuations in the SDW 
and SC channels are coupled at intermediate energies 
A > E > Ep between the bandwidth and the Fermi 
energy, but decouple at energies below Ep. The system 
behavior below Ep is governed by conventional ladder 
RG, and each vertex flows according to dTi/dL = Tf. 

For the toy 2-pocket model, earlier results showe d 22 i 51 
that SDW instability is the dominant one at perfect nest- 
ing. The SC vertex is repulsive at large energies but 
changes sign under parquet RG and become attractive 
above some RG scale. The SDW and SC couplings flow 
to the same value at the fixed point of RG equations, and 
the fixed point of parquet RG has extended 0(6) symme- 
try^ If the scale of Ep is reached before this fixed point, 
SDW order prevails at zero doping but is reduced and 
eventually destroyed at finite doping. Whether or not 
SC appears in place of SDW order depends on whether 
SC vertex already changes sign and becomes attractive 
at Ep. If superconductivity appears, the SC gap has a 
simple plus-minus structure. 

The main goal of this paper was to understand how this 
scenario is modified in realistic 4-pocket and 5-pocket 
models. We considered both models in the two limits: 
one when one of the two hole pockets centered at (0, 0) 
is weakly coupled to other pockets and can be neglected, 
and the other when the two hole pockets centered at (0, 0) 
are equivalent. We found identical results in both limit 
what gives us confidence that the system behavior in the 
intermediate case of two non-equivalent hole pockets at 
(0, 0) remains the same as in the two limits. 

Our main results are the following: 

• In both 4-pocket and 5-pocket models electron-hole 
and electron-electron interactions are generally an- 
gle dependent. The most relevant angle depen- 
dence is in the form cos 20, where <j) is the angle 
along an electron FS. Because of this angular de- 
pendence, there are three different vertices in the 
SC channel. One of these vertices turns out to be 
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attractive, in most cases, beginning from the largest 
energies. The symmetry of the attractive interac- 
tion is extended s±wave (as opposed to a conven- 
tional s + +). Other two SC vertices are repulsive 
at all scales. 

• This attractive SC vertex favors the s± state 
in which the gaps along hole FSs are angle- 
independent (up to cos 40 corrections), while the 
gaps along the two electron FSs are in the form 
A e ± A e cos 2(f). The interplay between A e and A e 
depends on the strength of cos 2</> component of the 
interaction and also on the interplay between intra- 
pocket and inter-pocket Coulomb repulsions. De- 
pending on the parameters, the electron gaps can 
be either nodeless (A e > A e ), or have accidental 
nodes (A e < A e ). 

• In 5-pocket model at perfect nesting, the SDW ver- 
tex remains larger than this attractive SC vertex. 
The two flow up to the same values at the fixed 
point, if this fixed point is at an energy larger than 
Ep, and the fixed point has enlarged symmetry. 
This behavior is exact when the vertices are angle- 
independent, but only very weakly changes due to 
angular dependence of the vertices. If the system 
flows down to Ep without yet reaching the fixed 
point, SDW order wins. Away from perfect nesting 
SDW order is suppressed, and the system eventu- 
ally develops a SC instability. 

• In 4-pocket model, the situation is similar at large 
E (i.e., at small RG parameter XogK/E), but before 
the fixed point of parquet RG is reached, SC vertex 
becomes larger than SDW vertex. If this happens 
before the scale of Ep is reached, the system devel- 
ops SC instability already at perfect nesting, and 
SDW instability does not appear. If SDW vertex 
remains the largest down to Ep , the system devel- 
ops SDW instability at and around perfect nesting, 
and SC instability at larger dopings. 

• We found that the SC gap is more likely to have 
accidental nodes on electron FSs in 4 pocket model 
than in 5-pocket model. Namely, for the same input 
parameters, there is a parameter range where the 
gap is nodal in 4-pocket model and no-nodal in 5- 
pocket model. This agrees with several RPA-type 



studies based on spin fluctuation o 20 ' 29 ' 30 . Still, we 
found that in both 4-pocket and 5-pocket model the 
gap can be either nodal or node-less, depending on 
parameters. 

Our analytical results are fully consistent with nu- 
merical fRG study of 4-pocket and 5-pocket models by 
Thomale et al 2 !. We view this agreement as the evidence 
that the differences between 4-pocket and 5-pocket mod- 
els are geometrical (different combinatorics in RG equa- 
tions), and are captured already within analytical one- 
loop RG. We note in this regard that we found that the 
difference between 4-pocket and 5-pockct models is not 
caused by the angular dependence of the interaction and 
holds even when interactions are angle- independent. 

The results for the 4-pocket case demonstrate that 
SDW order need not be pre-requisite to SC± order, al- 
though for most part of the phase diagram it does appear 
at perfect nesting, and SC only appears upon doping. We 
also emphasize that the interplay between SC and SDW 
is both " mutual support" and "competition". Namely, 
SC and SDW fluctuations tend to enhance each other, 
what is relevant is the fact that in the applicability range 
of parquet RG (when SC and SDW fluctuations talk to 
each other), both SC and SDW vertices diverge upon ap- 
proaching the fixed point. At the same time, SDW and 
SC orders compete with each othe r 58 ' 59 , meaning that SC 
order only emerges when SDW order is reduced enough 
by doping, and SDW order does not emerge at all if SC 
order emerges first already at perfect nesting. 
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